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ABSTRACT 



We present the analysis of RXTE monitoring data obtained during the January 2009 outburst of the hard X-ray transient lA 1 118-61. 
Using these observations the broadband (3.5-120 keV) spectrum of the source was measured for the first time ever. We have found that 
the broadband continuum spectrum of the source is similar to other accreting pulsars and is well described by several conventionally 
used phenomenological models. We have discovered that regardless of the applied continuum model, a prominent broad absorption 
feature at ~55keV is observed. We interpret this feature as a cyclotron resonance scattering feature (CRSF). The observed CRSF 
energy is one of the highest known and corresponds to a magnetic field of B~ 4.8 x 10^^ G in the scattering region. Our data also 
indicate the presence of an iron emission line presence that has not been previously reported for lA 1118-61. Timing properties of 
the source, including a strong spin-up, were found to be similar to those observed by CGROfBAISE during the previous outburst, but 
the broadband capabilities of RXTE reveal a more complicated energy dependency of the pulse-profile. 
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1. Introduction 

The hard X-ray transient lA 1118-61 was first disc overed dur- 

° ing a n outburst in 1974 by the ArielS satellite (lEyles et alJ 

] ll975l) . The outburst lasted for -10 days and no flux could 

. be observed after wards. Pul s ations with a period of 405.6 s 

■ were observed by lives et al] (1 19751) and were initially inter- 

I preted as th e orbital period of two compact objects. It was 

suggested by iFabian et aP (1 19751) that the observed period was 

due to a slow rotation of the neutron star. The optical coun- 

terpart was identified as the Be- star He 3-640/Wray 793 by 

IChevalier & Ilovaiskyl(ll975l) and classified as an 09.5IV-Ve star 



with strong B aimer e missio n lines and an extended envelope by 
iJanot-Pache co et all (1 19811). The di stance was estimated to be 



5 + 2 kpc ( Janot-Pacheco et aO 19811). The clas sification and dis- 
tance was confirmed by ICoe & Paynd (Il985l) by UV observa- 
tions of the source. The X-ray spectrum of the pulsar was fitted 
with a power law with a photon index of F ~ 1 , with a possi- 
ble spectral softening to F ~ 0.9 during the peak of the outburst 
(significant at Icr confidence level). 

A second ou tburst occurred in 1992 and was observed by 
CGRO/BATSE (ICoe et al.l 11994 . The measured peak pulsed 
flux was ~ 150 mCrab for the 20--100keV energy range, sim - 
ilar to the 1974 outburst (ICoe et al .111 9941: iMaraschi et al.lll976h . 
It was followed by a period of elevated emission ~ 2 5 days af- 
ter th e main outburst. This lasted for ~ 30 days (see ICoe et al.l 
11994 Fig.l). Pulsations with a period of 406.5 s were detected 
up to lOOkeV and the pulse profile showed a single, broad 
peak, asymmetric at lower (20-40 keV) energies. A spin-up 
of 0.016 s/day was observed during the decay of the outburst. 



The pulsed spectrum was described with a single-temperature 
optically-thin, thermal bremsstrahlung model with a temperature 
of (15.1 ± 0.5) keV for the main outburst and (18.5 ± 0.9) keV 
during the elevated emission. Multi- wavelength observations re- 
vealed a strong correlation bet ween the Hg, equ ivalent width and 
the X-ray flux, which allowed lCoe et al.l (11994 to conclude that 
expansion of the circumstellar disk of the optical companion was 
mainly responsible for the increased X-ray activity. The perias- 
tron passage would then trigger an outburst if enough matter had 
accumulated in the system. This conclusion was supported by 
the pulsations with a period of ~ 409 s, which were also detected 
in quiescence (iRutledge et al.ll200"7b . 

The source remained in quiescence until 2 009 January 4, 
when a third outburst was detected by Swift (iMangano et al.l 
l2009h . Pointed observations with Swift/XRT allowed the detec- 
tio n of pulsations w ith a period of 407.68 ± 0.02 s reported later 
by lManganol (l2009l) . The complete outburst was regularly mon- 
itored with RXTE. INTEGRAL observed the source after the 
main outbu rst and observed fla ring activity ~ 30 days after the 
main burst (iLevder et al.ll2009l) . Suzaku observed lA 1118-61 
twice, once during the peak of the outburst and also ~ 20 days 
later when the flux returned to its previous level (S. Suchy, in 
preparation). 

We report here on the monitoring observations by RXTE. A 
timing analysis to determine the pulse period of the source and 
the rate of the observed strong spin-up was carried out. Pulse 
phase averaged spectral analysis revealed an absorption feature 
at ~ 55 keV, confirmed independently with Suzaku (S. Suchy, in 
preparation). We interpret this feature as a cyclotron resonance 
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Fig. 1. Swift/BAT daily light-curve of the outburst in 
2009 January (dotted line). Observation times by various mis- 
sions are indicated with vertical lines. 



scattering feature (CRSF), which is observed for the first time in 
this source. 



2. Observations 

The outburst was observed by multiple X-ray missions. The 
source light-curve as observed by Swift/BAT with marked ob- 
servation times by various satellites is presented in Fig[T] The 
source flux peaked on 2009 January 14, (Swift/BAT countrate 
0.12 counts/s corresponding to ~ 500mCrab in the 15-50 keV 
energy range) and slowly decreased afterwards. The source 
was regularly monitored with RXTE between January 10 and 
2009 February 4, with a total exposure of 86 ks (PCA) and dead- 
time-corrected live time of 29 ks (HEXTE-B). The data from the 
HEXTE-A detector were not used as it was not rocking during 
the entire observation and the background could not be properly 
estimated. The data were reduced with the HEASOFT version 
6.8 and a set of calibration files version 20091202. The spec- 
tral modeling was performed with the XSPEC package version 
12.5. In. The energy range 3.5-25 keV was used for the PCA and 
1 8-120 keV for the HEXTE spectra. Based on an analysis of re- 
cent Crab observations performed during the same time frame 
(2008 December 16, - 2009 March 11, proposal ID P94802), a 
systematic of 0.5% was determined for the PCA data. No sys- 
tematic error was required for the HEXTE data. 



3. Timing analysis 

To determine the pulse period of lA 1118-61 the PCA 
light-curve covering the complete observation in the 3-50 keV 
energy range was used. We searched for the pulse pe- 
riod and pulse period derivative using the phase connec- 
tion technique (IStaubert et al.l 120091) . Our best-fit results give 
p,pin=407.719(9)s, Pspin = -4.6(2) x 10-^ss-\ at a folding 
epoch of MJD 54841.62 (the uncertainties given in parenthesis 
are at Icr confidence level and refer to the last digit given). The 
best-fit residuals for pulse arrival times are presented in Fig. O 
Note that the obtained parameters characterize the pulse pe- 
riod evolution during the RXTE observations only and are not 
consis tent with the pulse period value provided by Mangano 
(l2009l) . This inconsistence is anticipated, because the luminos- 
ity and hence the accretion rate were significantly higher dur- 
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Fig. 2. Best-fit residuals for the pulse arrival times determined 
using RXTE PCA lightcurve (i.e. diff'erence between pulse ar- 
rival times as observed (O) and as calculated (C) with an as- 
sumed pulse period and derivative). Solid, dashed and dotted 
lines correspond to fits including the pulse period derivatives up 
to second, first and none (constant period) respectively. 



ing the RXTE observations. The spin-up rate of the neutron star 
depends on the accretion rate, and therefore extrapolating the 
timing solution obtained close to the peak of the outburst (with 
a stronger spin-up) to the beginning of the outburst (where the 
spin-up rate is expected to be significantly lower) gives a longer 
pulse period than can be directly measured at the time. Note 
also that our solution does not account for the Doppler delays 
caused by orbital motion, because the parameters of the orbit are 
unknown and our data do not allow us to find an unambiguous 
solution for the orbit. The estimated value fo r an orbital pe riod 
from "pulse period"-"orbital period" diagram lCorbetl (Il986h lies 
in the 400-800 d range, which is much longer than the span of 
our data. The intrinsic spin-period evolution is expected to be 
complicated, so it is difficu lt to separate it fr om the effects of the 
orbital motion. Moreover, ICoe et al.l (1 19941) suggested that the 
orbit may be almost circular, and the outbursts may occur at any 
orbital phase, so a possible orbital period value derived from the 
comparison of outburst times is also potentially ambiguous. 

Both the spin-up rate ~ 0.04 sd"^ and ff ux are somewha t 
higher (by a factor of ~2-3) than reported by ICoe et"aD (1 19941) 
for the previous outburst observed with CGRO/BAISE. The 
stronger spin-up measured for the current outburst is likely 
caused by the higher ffux and consequently higher accretion rate 
(although it is difficult to directly compare BAT and BATSE 
ffuxes, because BATSE measured only pulsed ffux and the en- 
ergy ranges are slightly different). A potential difference in the 
orbital phase during the two observations may also be responsi- 
ble for the difference in the observed spin-up rate. 

A set of pulse proffles in several energy ranges was con- 
structed with the determined period. The pulse proffle signif- 
icantly changes with energy as shown in Fig. [3] At energies 
below lOkeV the pulse proffle has two peaks. The secondary 
peak amplitude decreases towards higher energies, disappear- 
ing above lOkeV. A shoulder appears on the other side of the 
main peak at about the same energy. The pulse proffle becomes 
single peaked and gradually more symmetric and narrow at 
higher energies. A similar behavior was observed previously 
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Fig. 3. Normalized "pulse phase"-"energy" matrix using PCA (below 22 keV) and HEXTE (above 22 keV) data and pulse profile 
evolution with energy. A slice at a constant energy gives a background- subtracted pulse-profile normalized to unity at the pulse 
maximum (shown in the middle pane, the pulse profiles are shifted with respect to each other to avoid confusion). Contours represent 
2cr and 3cr significance levels for the absorption feature in the residuals to the fit with CompTT model modified by photoelectric 
absorption and emission line at 6.4 keV. The pulse fraction, defined as (Amax - Afnin)/(Amax + ^mm) is shown as a function of energy 
in the right panel. 



with CGRO/BATSE (ICoe et al.|[T994l) . The pulse fraction, de- 
fined as (Afnax - Amin)l{Aniax + Amin) mcrcascs with energy in 
a similar way to other accreting pulsars, except for a drop at 
around 8 keV, where the second peak disappears (see Figl3]). 



4. Spectral analysis 

The longest observations were obtained during the maximum of 
the outburst, accordingly we focused on this data for the spec- 
tral analysis (observations 94032-04-02-03 to 94032-04-02-10). 
The spectra of lA 1118-61 in diff'erent energy ranges were 
previously described by power law an d bremsstrahlung models 
(iMaraschi eraDll976l:[Coe et al.lll994 . but we found that these 
models do not describe our data adequately. Our results show 
that the broadband continuum of lA 1118-61 ca n be well de- 
scribe d by the FDCU T, NPEX, and CompTT m odels (ICoburn et al.1 
120021: lMiharalll995l:lTitarchuk & Lvubarskiilll995h with qualita- 
tively similar goodness of the fit, modified by photoelectric ab- 
sorption and iron emission line at ~ 6.4 keV. 

All used continuum models in their best-fit residuals show a 
prominent absorption feature at ~ 55 keV. The inclusion of an 
absorption line with a Gaussian optical depth profile leads to a 
significant reduction of ;^^ in all models (see Table [T]). We in- 
terpret this absorption feature as a CRSF, which is detected in 
the spectrum of 1 A 1 1 18-61 for the first time. The observed en- 
ergy of the feature is one of the highes t known and correspo nds 
to a magnetic field B~ 4.8 x lO^^G (iHarding & Laill2006h in 
the scattering region. Preliminary phase-resolved analysis shows 
that the energy of the CRSF is not likely to change significantly 
with the pulse phase (see FigO. 

The best-fit residuals of FDCUT and NPEX models show 
also a prominent emission line-like feature at ~ 8 keV. Similar 
features at diff^erent energies were reported f or a number of 
sources (ICoburn et al.l 12002 : iRodes-Roca et al.l [2Q09 and refer- 
ences therein). ICoburn et aP (l2002h suggested that the employed 



phenomenological models for the continuum may be oversim- 
plified for the real source s and hence may be resp onsible for this 
eff'ect. On the other hand. lRothschild et all (l2006l) suggested that 
the feature may be associated with the fluorescence copper line 
from the Be/Cu collimator of the PCA instrument. To clarify the 
situation a more detailed analysis of the Suzaku data is currently 
made and will be presented as a separate paper (S. Suchy et al, 
in preparation). 

We found that the inclusion of an emission line with 
Gaussian profile may help to account for this feature and does 
not significantly aff'ect other model parameters. The;^^ substan- 
tially improves for fits with FDCUT and NPEX models, accord- 
ingly we included the line in all fits with those models. The 
energy and the width of the line were fixed at 8.04 keV and 
0.01 keV, correspondin g to the copper K^ line, as proposed by 
iRothschild 



orresp onamg 
etal.f(l2006h . 



The fit may be slightly improved by the inclusion of an ad- 
ditional absorption feature at ~ 23 keV for PCA data only. The 
investigation of the Crab residuals (proposal P94802) shows that 
a similar feature is also present in the PCA Crab spectrum (more 
than 99% s ignificance with max imum likelihood ratio (MLR) 
and F tests (iProtassov et al.ll2002l p-value ~ 0.006). This is con- 
sistent with our lA 1118-61 data, because the line is not re- 
quired by HEXTE data. On the other hand, a shallow funda- 
mental line at 22keV (with a first harmonic at 55 keV) could 
be more difficult to detect due to this instrumental feature using 
PCA data. A comparison of our RXTE results with the Suzaku 
observations will help to clarify the picture (S. Suchy et al, in 
preparation). 



5. Conclusions 

For the first time since 1992 a major outburst of the Be/X- 
ray binary lA 1118-61 was observed in 2009 January with 
RXTE, following the trigger of Swift/BAT. Strong pulsations 
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Fig. 4. Best-fit residuals for various continuum models of lA 1 118-61. A line-like emission feature at ~ 8 keV was modeled with a 
narrow line with Gaussian profile as described in text. The absorption feature at ~23 keV was included for PCA only, to account for 
the similar residuals in the Crab spectrum (see text). 
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Table 1. Best-fit results for diff'erent models. An emission line {Cu) with energy and width fixed at 8.04 keV and 0.01 keV was 
added to FDCUT and NPEX models. All models include also an absorption like feature at ~ 23 keV for PCA to account for similar 
residuals seen in Crab spectra (with the line energy and width fixed to those obtained from Crab fits). For CompTT model the;^^ may 
be improved by allowing the line parameters to vary (last column). 



Notes. (^) Parameter frozen during the fit. ^^^ [keV] ^'^ [lO'^ phkeV" 



1] (^^ [10-3 phcm-2 s-i] ^'^ [atoms cm-2] 



with a period of 407.72 s and a spin up of -4.6 X 10"^ SS"^ Acknowledgements. VD thanks DLR for financial support (grant 50OR0702). 
were detected. A similar temporal behavior was observed with 
CGi^O/B ATSE during the previous outburst. A broadband spec- 
trum of the source was obtained for a first time, and an absorp- 
tion feature at ~ 55 keV, interpreted as a CRSF was detected. The 
inclusion of the feature significantly improves fit results with all 
applied continuum models and its energy does not depend sig- 
nificantly on the model used. More detailed spectral analysis is 
ongoing and the results will be published elsewhere. 
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